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We conduct numerical simulations in Painleve-Gullstrand coordinates of a variety of K-essence- 
type scalar fields under spherically symmetric gravitational collapse. We write down generic con- 
ditions on the K-essence lagrangian that can be used to determine whether superluminality and 
Cauchy breakdown are possible. Consistent with these conditions, for specific choices of K-essence- 
type fields we verify the presence of superluminality during collapse, while for other type we do not. 
We also demonstrate that certain choices of K-essence scalar fields present issues under gravitational 
collapse in Painleve-Gullstrand coordinates, such as a breakdown of the Cauchy problem. 

I. INTRODUCTION 

In recent years, there have been a number of models put forward to explain the accelerating expansion of the 
universe. Falling under the general rubric of dark energy, such models have included quintessence [1 , Elko spinors 
[2], H-essence [3], quintom matter [4], phantom matter [5] and more. However, these models do not explain why the 
expansion of the universe is accelerating while humans are here to observe it (the cosmic coincidence problem). Models 
that avoid anthropic rationales for explaining cosmic coincidence generally require the fine-tuning of parameters. It 
was with the intention of avoiding both fine-tuning and anthropic arguments that K-essence [SJ [7] was first proposed. 

K-essence models invoke dark energy in the form of a scalar field with non-linear kinetic energy. These models 
are constructed so that the scalar field develops a negative pressure once the matter dominated era begins, and so 
the associated dynamical behaviour is then deemed responsible for the accelerating cosmic expansion of our universe. 
K-essence initially came under criticism [8] because it was shown that any such model that solved the coincidence 
problem also necessarily involved perturbations propagating faster than light. However, the proposers of K-essence 
defended their model [9] by demonstrating that despite superluminality, causality was preserved for physical situations. 

It is due to this intrinsic superluminality that the idea of gravitational collapse of K-essence matter has become 
recently popularized. A number of K-essence-type models, characterized by non-linear kinetic energy, have emerged. 
These models serve the primary purpose not of solving the coincidence problem, but of exploring the idea of super- 
luminal information escaping from a black hole. It was shown initially in |10j that for the stationary black hole case, 
signals could escape from within the light horizon as superluminal perturbations. More recently, a dynamic collapse 
scenario was considered and numerically simulated in maximal slicing coordinates |llj . This study found evidence of 
superluminality in certain Lagrangian choices for K-essence. Other choices of Lagrangian, however, showed no such 
evidence. This study also demonstrated that in some choices of Lagrangian we are limited by the eventual breakdown 
of the Cauchy problem. 

Because the well-posedness of the Cauchy problem is dependent on the time-coordinate remaining globally valid, the 
authors of [TT] were unable to determine whether the observed breakdown was a result of the choice of Lagrangian or 
of coordinate systems. To this end, we now examine a selection of choices of Lagrangians in Painleve-Gullstrand (P-G) 
coordinates, which are regular across future luminal horizons. We provide general analytic expressions for determining 
whether superliminal propagation and/or Cauchy breakdown occur for specific classes of K-essence fields. Because 
P-G coordinates offer an indisputably distinct time coordinate from maximal slicing coordinates, simulations in P-G 
should illuminate the coordinate dependence of Cauchy breakdown or lack thereof. Using P-G coordinates is also a 
strong check to ensure the coordinate independence of the results of the maximal slicing simulations in general. 

Through a series of numerical simulations in P-G coordinates of three models of K-essence, we demonstrate that 
the gravitational collapse of K-essence can produce either superluminal or subluminal sonic horizons. Hence we verify 
the qualitative results of [TT] and therefore provide evidence for the coordinate independence of these results. We 
further examine the breakdown (or lack thereof) of the Cauchy problem for these three models. We find that one 
model experiences a breakdown of this problem in many parameter choices, and that the other two models do not 
experience Cauchy breakdown. This corroborates the results found in [TT]. Therefore, our results are indicative that 
the issue of Cauchy breakdown for these models of K-essence is coordinate independent. 
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II. METHODS 



We start with a general action for K-essence coupled to gravity, given by 

7 = /^ 4 Kt^ + £P0 ) (1) 

where is the spacetime scalar curvature, G is Newton's gravitational constant in four dimensions, and X = 
— ^W a ipW a tp where ip is the K-essence scalar field. 

If we vary the action with respect to the metric, we obtain the Einstein field equation 

Gab — KT ao (2) 

where we have denoted denoted the Einstein tensor as G ao , and we note that for K-essence we have 

Tab = Cx^aip^bip + Cg ab (3) 

Above and throughout this paper we use Cx to mean dC/dX and Cxx to mean d 2 C/dX 2 . We vary the action with 
respect to ip, the K-essence scalar field, and obtain equations of motion for K-essence: 

.g a "V a V^ = (4) 

Note that in the above we have used g ab , which is the effective (inverse) metric governing propagations of perturbations 
of the K-essence field, given by: 



g ab = C x g ab - C xx V a 1>V b 1> =4> g ab = ^g ab + c^%^V a ^V 6 ^ ( ) 

Cx L x 

where 



2 £x 



C x + 2XC X x 
is the speed of sound. 

The luminal apparent horizon forms when surfaces of constant r become null, implying 

g ab V a rV b r = 0=^g rr =0 (6) 

The right hand side assumes only a coordinate basis with r as the spatial coordinate. 

We are now able to examine general conditions for superluminal propagation and Cauchy breakdown. The condition 
for the formation of the sonic horizon requires that the same general condition as ^ be satisfied for the effective 
metric of K-essence, that is 

g ab V a rV b r = (7) 

which, in analogy to the luminal case leads to: 

9 rr - (8) 

A calculation reveals that in general 

§rr = Cxfr _ c ^ a = + 2XCxx) _ Cxxj^-Aii) 2 = - Cxxj-^Aii? = (9) 



where in the last line g is the determinant of the 2- metric as given in equation (24 1, and I is an arbitrary length scale. 
Note that g is negative, so — g is positive. Eq.Q shows that the relative sign difference between the sonic horizon 
condition and the luminal horizon condition is determined by the signs of Cx and Cxx ■ Superluminal escape requires 
the sign of g rr to be positive just inside a black hole horizon, where g rr is negative. We assume that Cx is positive, 
which is necessary for the kinetic term in the stress energy tensor to be positive (cf. Eq.Q). Thus, for the right hand 
side of equation ([£]) to be positive, it is necessary (but not sufficient) that Cxx < 0. When this condition is satisfied, 
we can generically say that superluminal escape is possible. This will be verified in the numerical studies below. 
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Another key aspect of our method is the examination of the potential breakdown of the Cauchy problem. A 
breakdown of the Cauchy problem occurs when the surfaces of constant time ^ become null with respect to the 
effective metric. We can write this as 



g ab W a W b t = => g u = (10) 



which using ( 39 ) , may be suggestively rewritten as 



g" = C x g tt - L xx (W) 2 = ^fg u - ^x-^-^'f (11) 

with g and I defined as in equation pi). In simulating equation ([2]), we require that surfaces of constant time be 
spacelike; in other words we require that g < 0. However, for K-essence we must then require that in addition to 
this, the surfaces of constant time of the effective metric are also spacelike. That is, g u < 0. Clearly, for certain 
Lagrangian choices, g u can go to zero independently of g tt . When surfaces of constant time are no longer spacelike 
with respect to either the background or effective metric, we say that the Cauchy problem ceases to be well posed 
on this surface, and our simulations must halt. Eq.(ll) and Eq.([9| show that the condition for Cauchy breakdown 



in a theory is the same as those for superluminal propagation, though the onset of when each occurs depends on the 
relative magnitude of \g rr / g lt \- 

Of course the above assumes that both C x and C xx are everywhere regular. It is also possible for pathologies to 
occur if either becomes singular. This too will be illustrated in one of the examples to follow. 

In order to be able to distinguish physical properties of the K-essence field from a breakdown of the coordinates it 
is important to use spatial slicings that are regular at luminal horizons. To this end we work in Painleve-Gullstrand 
coordinates. Our simulations are a generalization to the K-essence case of the method developed in [T2]. In Painleve- 
Gullstrand coordinates the metric take the form: 



ds 2 = -a 2 dt 2 + f dr + y ^^adt) + r 2 dVL 2 (12) 

where in the dynamical case both the lapse function a and the Misner-Sharpe mass function M. are functions of both 
r and t. Note that spatial slices (dt = 0) are flat so that P-G coordinates are singular only at the origin, r = as 
long as a 2 and M. are positive. This means that we are able to continue to evolve the simulation inside any horizons 
that form (as opposed to the situation in Schwarzschild coordinates, for example). Another reason for our choice of 
coordinates is simply to avoid issues associated with the usual choice of polar-radial coordinates. This choice would be 
inappropriate in this setting, because such coordinates are only valid until the first trapped surface forms. In this case, 
since we have superluminality potentially present, K-essence is not trapped by normal trapped surfaces, but rather 
by the trapped surfaces of the effective metric. By working in P-G coordinates we instead can follow the evolution 
of the system much further, providing a different qualitative perspective from that of maximal slicing coordinates, as 
noted above. 

Since we are considering spherically symmetric collapse, the number of dynamical degrees of freedom is reduced to 
that of tp and its conjugate momentum 11.0. These quantities obey the equations ( 39|40 ) given in appendix A, which 



contains more detailed description of our dynamical variables and specific numerical methods. For now, it will suffice 
to note our initial conditions on ip and II,/,, which we take to be 

4(t = 0) = Ar 2 exp ( (13) 

as illustrated in fig [T] and 

n^(* = o) = o (u) 

for n^,. For all of our simulations, we chose B = 0.3 and r = 1.0. 



III. RESULTS 



Before proceeding to the results for each specific choice of Lagrangian, we review the behaviour of a spherically 
symmetric gravitational collapse scenario for 'normal matter' (where C(X) = X) and note the conditions for the 



FIG. 1: V vs r at t = 



formation of light apparent horizons and sonic horizons (horizons trapping the perturbations of the effective metric) 
in P-G coordinates. 

Let us recall that for initial data such as that we have chosen, there are a number of different possibilities for how 



the collapse of normal matter proceeds. If the amplitude of the scalar field, A, as given by equation ( 13 ) is sufficiently 
small, subcritical collapse occurs. This means that no horizons are formed, and the ingoing pulse becomes, at the 
origin, an outgoing pulse and the field simply 'passes through' itself. However, if A exceeds some critical value, the 
collapse becomes supercritical. Luminal apparent horizons form, from inside which nothing can escape. 

K-essence matter is governed by the effective metric g, and as shown above perturbations in this matter are 
permitted to exceed the speed of light. In P-G coordinates, the two distinct conditions for the formation of the 
luminal and sonic horizons, are respectively 



r - 2GM = (15) 



and 




r J (C X +2XC XX )~C XX \ 1 ^J Z - + X I— 1 P') = (!(,, 



Similarly we can find a general condition for the breakdown of the Cauchy problem in P-G coordinates. By noting 
that this occurs when surfaces of constant t are no longer timelike (as stated above), we obtain the general condition 



( n v 2 



^+( 7 -J F -) C XX =Q (17) 



Eqs.( 15|l6|17 1 are used in the numerical code to locate luminal horizons, sonic horizons and Cauchy breakdown, 



respectively. 



A. Born-Infeld Lagrangians 

We examine two choices of Born-Infeld Lagrangians. These are based on the Lagrangian given in [TU], which 
was proposed to demonstrate information escaping from inside a stationary luminal apparent horizon. As such they 
represent interesting models for studying gravitational collapse since it is not immediately clear what effect the 
superluminal character of collapsing K-essence matter will have on this process. 

The first of this type of Lagrangian is given by 



£ 1 (X) = a 2 [ ] Jl + — -1]-A (18) 

where a is the K-essence parameter. As a goes to oo, we recover C(X) = X, the standard Lagrangian for a massless 
scalar field ('normal' matter). We choose A = for our simulations. Note that for gravitational collapse it is possible 
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FIG. 2: Sonic (dashed) and light (solid) horizon formation conditions vs r at P-G time t = 0.172499478 for &{X) with A=0.13 
and a — 3.0 



for X to become negative, unlike cosmological scenarios where homogeneity implies all physical quantities depend 
solely on time ensuring X > 0. 

From figure [2] we see that for the Lagrangian given in (18), we have the possibility of superluminal escape from the 



black hole. For the given choice of parameters, we clearly see that the sonic horizon forms inside the apparent horizon, 
indicating that perturbations in the intermediate region which travel faster than light can escape the light horizon. 
The spatial distinction between the two horizons does not persist for long though: by P-G time t = 0.199998676, 
the two horizons have merged spatially within the precision of our numerics. This behaviour was also qualitatively 
observed in maximal slicing coordinates in |llj . As described therein, the reason for the rapid merger of the sonic 
and apparent horizons is that the K-essence scalar field quickly falls into the black hole or propagates outwards away 
from the horizon. This leaves the background and effective metrics equivalent in the horizon region, so luminal and 
sonic horizons are both dependent only on the background metric. This behaviour is equivalently demonstrated in 
figure [3J which depicts the null geodesies of both the background and effective metrics, as well as both the sonic and 
luminal apparent horizons for a typical choice of parameters. We find, somewhat surprisingly, that the inner sonic 
and luminal trapping horizons remain separated, while the outer ones merge very rapidly, leaving only a very narrow 
region of superluminal escape. 

Aside from the striking example provided in figure [2] it will be instructive to provide a more in-depth discussion 
of the behaviour of this model in the parameter space explored. This will also serve as a discussion of breakdown of 
the Cauchy problem in this model. We examined this choice of C for a range of parameters: our values of A ran from 
0.01 to 0.15, and our values of a from 0.05 to 400. Within this range, there were a number of different categories 
of behaviour observed. In general terms, all parameter pairing for this choice of C may be grouped into one of the 
following categories: 

1. Subcritical Collapse - no Cauchy breakdown 

2. Cauchy breakdown prior to any horizons forming 

3. Cauchy breakdown after apparent horizon formation but before sonic horizon formation 

4. Cauchy breakdown after both horizons form, where sonic horizon appears inside light horizon 

5. Cauchy breakdown after both horizons form, where within our numerical precision, sonic and apparent horizons 
are indistinguishable 

Before discussing specifically for which parameters each behaviour is present, we note the factors which could create 
different behaviours in different parts of the parameter space. First, we know that as a increases we approach the 
behaviour of 'normal' matter, and as a decreases we depart from this behaviour. Therefore, we expect breakdown of 
the Cauchy problem to occur sooner for cases of smaller a. Secondly, we know that as A increases, the collapse will 
go from from subcritical to supercritical. Within the supercritical range, as A becomes bigger, collapse will be more 
rapid and horizon formation will occur sooner. With this in mind, we summarize our observations for C\ as follows. 



1. Subcritical collapse without Cauchy breakdown occured for A = 0.01 — 0.04 for relatively large a - the minimal 
value being between 0.01 and 2.0, depending on the value of A, and extending to a = 400.0 in all cases. This is 
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FIG. 3: Map of null geodesies for Ci(X) of the background metric (dark blue solid) and the effective k-essence metric (dark 
green dashed) as well as the luminal apparent horizon (light blue solid) and sonic apparent horizon (light green dashed), for 
parameters A=0.12, a=2.0. The red line is the sonic event horizon, i.e. the last outgoing sound wave that does not quite escape 
the sonic trapping region. The luminal event horizon is just slightly to the right of the sonic event horizon, but is too close to 
distinguish on this scale. 



as we expect 
collapse. 



large a implies normal matter behaviour, and small A in normal matter corresponds to subcritical 



2. The breakdown of the Cauchy problem prior to any horizon formation occured at small a for all values of A. 
Figure [4] shows an example of this type of behaviour. Depending on the value of A, the maximum value of a for 
which this occurred was between 0.075 and 80.0. This is explained by the fact that, as stated before, at small 
a we have behaviour less typical of normal matter and we expect Cauchy breakdown to occur sooner. 

3. We observed 2 cases of Cauchy breakdown occurring between the formation of the sonic horizon and the forma- 
tion of the luminal apparent horizon. This was for A — 0.11 and a — 0.5,0.75. This behaviour is understood 
easily once we understand that in this choice of C sonic perturbations can be superluminal - the luminal apparent 
horizon would then of course form before the sonic apparent horizon. In these cases, Cauchy breakdown occurs 
between these formations. 

4. For intermediate values of a and large values of A, we observe breakdown of the Cauchy problem after luminal 
and sonic horizons have formed, with the sonic horizon initially occurring inside the luminal horizon. This type 
of collapse was observed for A = 0.1 — 0.13 and for values of a ranging from 0.5 to 25.0 depending on the value of 
A. As we increase a within this range, the window for superluminal escape becomes smaller in time, or in other 
words the two horizons converge more rapidly. This is as expected, since as a increases we progress towards 
normal matter behaviour, so Cauchy breakdown occurs later and the two horizons become indistinguishable. 
The two horizons also converge more rapidly as A is increased, because increasing A increases the rapidity of 
the collapse. 

5. For large values of a and large values of A, we find that Cauchy breakdown occurs after both horizons have 
formed. However we are unable to numerically distinguish the locations of the horizons to within our levels of 
tolerance. Again, this is as expected because for normal scalar field collapse there is no distinction between the 
luminal and sonic horizons and large a is the normal matter limit. 

Summarizing, for this choice of C we observe that superluminal escape is possible for a portion of the parameter 
range. To observe this, we must have a sufficiently large to allow us to observe horizons prior to Cauchy breakdown, 
and we must also have a sufficiently small so that that the two horizons are numerically distinguishable. We also 



FIG. 4: Cauchy breakdown condition vs r for Ci at P-G time t=0. 497767983 (the last timestep before Cauchy breakdown) 
for A=0.1 and a = 0.5. The point of closest approach to zero is approximately r=0.3434. 




FIG. 5: Sonic (dashed) and light (solid) horizon formation conditions vs r at P-G time t = 0.744807261 for C%{X) with A=0.1 
and a = 2.0 



require A large enough to avoid subcritical collapse, but small enough to avoid the two horizons converging too rapidly 
to be distinguished numerically. 

We also see for this choice of C that we have Cauchy breakdown occuring at small values of a for all A, in accord 
with results for the MS case [TT] . That we observe this breakdown for similar regions of the parameter space using an 
entirely different time coordinate indicates that the breakdown of the Cauchy problem in this choice of Lagrangian 
is likely a generic coordinate-independent feature. Conversely, for A = 0.05 — 0.13, the breakdown of the Cauchy 
problem occurred eventually for all values of a which we examined. Therefore we see that in P-G coordinates the 
breakdown of the Cauchy problem is not restricted to small a as it is in maximal slicing coordinates. 



We can easily modify the Lagrangian ( 18 ) to avoid the problem of Cauchy breakdown by choosing a second Born- 
Infeld type Lagrangian [TT] 



C.iX)=n-[l~ \ l-^]-A (19) 



where again taking the limit as a goes to oo yields normal matter. 

Setting A to 0, in this case we find no evidence in our simulations of superluminal behaviour. Our range of 
parameters observed in this case was the same as for L\. For a typical choice of parameters we see in figure [5] that 
the sonic horizon forms outside of the light horizon, as we would expect for normal matter. Hence, there is no 
superluminal propagation in this case. This is further illustrated in figure [6] by mapping null geodesies and horizons 
for both background and effective metrics. As for Ci, we observed the two outer trapping horizons merging quickly, 
for the same reason, and the inner ones remaining separated. For the parameter values used in figure [5j the two 
horizons were indistinguishable by P-G time t=0. 779575347. 




FIG. 6: Map of null geodesies £.2(X) of the background metric (dark blue solid) and the effective k-essence metric (dark 
green dashed) as well as the luminal apparent horizon (light blue solid) and sonic apparent horizon (light green dashed), for 
parameters A=0.1, a— 2.0 



While we did not see any Cauchy breakdown in this case, we did encounter another issue with this choice. As can be 



clearly seen by the form of equation (19), it is possible for Cx = 



to become infinite when X 



However, 



[9] clearly states that Cx should not be allowed to become infinite. Our simulations indeed become numerically 
unstable in cases where X evolves to this value. 

In the MS case [TT], simulations were reported to eventually fail not for this reason, but due to the formation of 
shock waves of perturbations. This is indicated by the values of the speed of sound c s going to 0, where 

2 £x 



C 



x 



2XC 



1 



xx 



2X 



(20) 



and so c s — > as X approaches \. 

The full description of behaviour for each portion of the parameter space we examined is analogous to that for C\ 
in many ways, so we present a slightly abbreviated version of our results compared to C%. 

1. Subcritical collapse: as before, this occurred for small A and relatively large a, consistent with the limit of 
normal matter. 

2. We found that Cx could become infinite after the formation of only a sonic horizon. This is somewhat analogous 
to case 3 for C\, in which Cauchy breakdown occurred after the formation of only a luminal horizon. In this 
case, Cx becomes infinite at small values of a for all values of A. We note that in reality this case encompasses 
two distinct cases. In one case, we have this occuring for values of A which at larger a produce both horizons 
prior to Cx becoming infinite (see below). This case is not unexpected: we simply see Cx diverging prior to the 
expected luminal horizon forming. The second case is more thought provoking: this behaviour also occurs for 
values of A for which at larger a we see subcritical collapse with no horizons of any type forming (see above) . 
This indicates there must be a threshold value of a below which subcritical behaviour yields to the formation 
of a sonic horizon. This behaviour is unexpected and may merit further investigation. 

3. For mid-range values of a and large values of A we found that Cx became infinite after both horizons formed, 
with the sonic horizon appearing outside the light horizon. As before the two horizons were visibly separated 
initially, then converged together. 

4. Finally, for large a and large A, we again found both horizons forming prior to Cx becoming infinite. We were 
unable to distinguish between the two horizons within our numerical limits of tolerance. 
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FIG. 7: Sonic (dashed) and light (solid) horizon formation conditions vs r at P-G time t = 0.399778586 for £3 (A") with A=0.12 
and C=0.8 



As before, we observed larger A to imply more rapid gravitational collapse and hence horizon formation taking place 
sooner. While we found, as expected, no Cauchy breakdown in the parameter space for this model, we encountered 
a new problem: that of diverging Cx ■ Specifically, we found that Cx became infinite more rapidly for smaller values 
of a, consistent with increasing departure from the normal matter limit. Our results are commensurate with the MS 
case [TT] and indicate that the Cauchy problem (or lack thereof) and the problem of diverging Cx (or vanishing c s ) 
are likely coordinate independent. 



B. Lagrangians Allowing a Global Time Coordinate 

We next examine the behaviour of a choice of Lagrangian that is completely free from the problem of Cauchy 



breakdown by construction. Such a Lagrangian would have the left hand side of (17) be either positive everywhere 
or negative eveywhere. One such example is 

C 3 {X) = ^L+e&-l (21) 
where C is a parameter of the model; as C goes to oo we recover the behaviour of normal matter (C = X). We can 



verify quite easily using a more general form of (17) that this makes the left hand side of (17) positive everywhere, 
therefore entirely avoiding Cauchy breakdown. 

Upon simulation of this model, we do not find evidence of superluminality. Commensurate with the MS case 
we observe in certain parameter choices the sonic horizon forming outside of the apparent horizon, as in the case 
of normal matter and for C^(X). We illustrate such results for a typical choice of parameters in figure [7| This 
offers an indication that the subluminality of this particular global-t model is independent of the coordinate choice. 
We also see the same behaviour as in the previous two cases with respect to the eventual merging of the horizons. 
For the parameter values chosen in figure [7J for example, we see the two horizons as indistinguishable by P-G time 
t=0. 449391824. This is perhaps best illustrated in the map of null geodesies and horizons given in figure [8] We do 
not find Cauchy breakdown in this expected by construction. 

For this model, we examined parameter pairings in the range of A = 0.02 — 0.13 and C = 0.5 — 100.0. We note that, 
as discussed in [11], the hyperbolicity requirement of equation Q prevents us from choosing C < 0.45. Although we 
encountered some numerical difficulties (explained below) that forced the program to halt before completion, these 
were unrelated to the physical or mathematical characteristics of the model, in contrast to the situation for C\ and 
£2- With this in mind, we give the following overview of different types of behaviour found for this case. 

1. Subcritical collapse was observed for small A for all choices of C. 

2. For mid-range values of A and small values of C, we observed a sonic horizon forming prior to the code halting, 
but not a luminal horizon. This is as expected from the construction of the model. 



3. For mid-range values of A and larger values of C, we observed neither sonic nor luminal horizons forming before 
the code halted. This is due not to a different run time before halting (which was comparable at different values 
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FIG. 8: Map of null geodesies £z(X) of the background metric (dark blue solid) and the effective k-essence metric (dark 
green dashed) as well as the luminal apparent horizon (light blue solid) and sonic apparent horizon (light green dashed), for 
parameters A=0.12, C=0.8 

of C) but rather due to the fact that at higher C (nearer the normal matter limit) the two horizon formations 
occur more closely together in time, so neither form before the code halts. 

4. For large values of A and small values of C, we observed subluminal behaviour such as that shown in figure 
[7j with the sonic horizon forming outside the luminal horizon. This is a natural extension of the behaviour at 
midsize A and small C. Here, since A is larger, the horizons form earlier, and so both have time to form before 
the code halts. 

5. For large values of A and large values of C, we see both horizons forming before the code halts, but they are 
indistinguishable within our degree of precision. 



We pause to comment on the numerical difficulties we encountered for this case. Because of the recursive nature 
of the definition of X in P-G coordinates (see appendix A), we were forced to use a numerical solver to ascertain 
the value AT at a given time-step. We chose a Newton-Raphson solution method, which had a very good rate of 
convergence. However the tendency of X to become unpredictably large near the origin after horizon formation made 
it extremely difficult to find an initial guess that would yield convergence for the entire duration and spatial extent 
of the simulation. Hence we were forced to halt the code on those occasions that the Newton-Raphson method failed 
to converge. 



IV. CONCLUSIONS 



Through numerical simulations of the gravitational collapse of K-essence scalar fields in Painleve-Gullstrand coordi- 
nates, we have been able to observe a number of interesting phenomena. For the Lagrangian C\ we have numerically 
verified that sonic horizons form inside the luminal apparent horizon, allowing superluminal perturbations to escape 
the black hole. For other choices of Lagrangian (£2 and £3), we see that collapse proceeds with the sonic horizon 
outside of the luminal horizon, indicating no possibility of superluminal escape. For all three Lagrangians we studied 
the two different horizons converge relatively quickly, as the K-essence scalar field vacates the region of the horizons. 
These results are fully commensurate with the MS case [TT], suggesting that they are quite general and coordinate 
independent. 

We can understand these features by examining the general form of the sonic horizon condition (16). Superluminal 
behaviour generically happens when the luminal horizon forms before the sonic horizon, yielding a region in which 
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(at least initially) light is trapped but not sound. As explained above, when Cx is positive this is definitely possible 
only when Cxx is negative in the trapping region. For £1, we have: 

£hx = 7- = = f ^f £i,xx = 9/1 , „„, oxo/ 9 (22) 



V/1 + 2A7«2 a 2 (1 + 2X/a 2 )3 /2 



whereas for 



^2,x — / =fff £2, xx — n „, (23) 

y/l - 2X/a 2 a 2 (l-2X/a 2 ) 3 / 2 

Moreover, for £3 all three quantities are manifestly positive. Thus, superluminal escape is in principle only possible 
for £1, and this is born out by our numerical calculations. 

It is important to note that the above considerations are slicing independent. The presence, or not, of superluminal 
propagation inside a trapped surface of the luminal metric is coordinate invariant, although the precise location of 
the trapping horizons may depend on the slicing. 

The possibility of superluminal escape from a black hole can be considered in terms of the sonic and luminal event 
horizons of the black hole. These are defined as the last sound and light rays, respectively to just barely escape from 
the black hole. The sonic and luminal event horizons coincide along a line of constant r, say r^, in the future after 
all the matter has fallen through. For C\ the sonic trapping horizon is to the future of the luminal trapping horizon, 
as shown in Figure [3j In this case, as one moves backward in time the sonic event horizon (red dashed line) "peels 
off" towards small r slightly faster than the luminal event horizon. Thus, there will be some K-essence perturbations 
just outside the sonic event horizon that will be able to escape to infinity from within the luminal event horizon. A 
related question is whether or not sonic perturbations are able to escape from inside the luminal trapping horizon 
(light blue solid line) . Although it is impossible to distinguish on the scale of Figure [3j we have verified numerically 
that as the sonic event horizon separates from the trapping horizon, it does enter the luminal trapping horizon. For 
the parameters that we have considered the corresponding region of superluminal escape is very small. It is not clear 
whether or not there is a range of initial data that allows a region of superluminal escape that is large enough to 
affect the information loss problem. This is a question that is worthy of further investigation. 

Note that none of these three choices of Lagrangian are suitable for addressing the cosmic coincidence problem [3] . 
For C\ we find that X becomes negative in certain regions of space, in contradiction with the assumed homogeneity 
of cosmology [11] . In the cases of £2 and £3, superluminality is not observed in either the MS or P-G simulations and 
so these models cannot be used to create the dynamical behaviour used to solve the problem of cosmic coincidence 
[5]. Any choice of £ that solves this problem must incorporate superluminality. These models are 'K-essence' only in 
the sense of having non-linear kinetic energy. 

As long as Cx and Cxx are regular, only C\ permits a breakdown of the Cauchy problem and this too was verified 
by the numerical calculations. However, in the case of £2, when X is too large, Cx and Cxx blow up. Our simulations 
(and those of [IT]) verified that this behaviour does indeed occur in spherically symmetric collapse with £ 2 . 

For £1, we have observed breakdown of the Cauchy problem. The smaller the value of a, the earlier this breakdown 
occurs during the simulation. This indicates that the constant time surfaces in P-G coordinates fail to continue to be 
Cauchy surfaces for either the background metric or the effective metric of K-essence, a phenomenon also observed for 
the maximal slicing time coordinate |11) . This provides a compelling (albeit inconclusive) indication of the coordinate 
independence of this generic breakdown of the Cauchy problem for these choices of £. Whether or not avoidance 
of Cauchy breakdown is necessarily connected to subluminality remains an open question; it remains to be seen 
whether another time coordinate may offer a well-posed Cauchy problem throughout for such choices of C. If no time 
coordinate exists that ensures the Cauchy problem is well-posed throughout collapse for both the background and 
effective metric, this would invalidate the predictive power of such models |10j . and eliminate the possibility of the 
physically sensible escape of information from a black hole. 



Acknowledgements 

We thank David Garfinkle and Ryo Saotome for helpful discussions and correspondence. We are grateful to 
SHARCNET for providing computer resources. This work was supported by the Natural Sciences and Engineering 
Research Council of Canada. 



12 



Appendix A 

In the following, where there is potential for confusion we denote 2 dimensional quantities with a hat and leave 4 
dimensional quantities as they are. Since we are considering spherically symmetric collapse, we dimensionally reduce 
the metric from 4 dimensions to 2 via 

ds 2 = -j—g^dx^dx" + r 2 (dO 2 + sm 2 (e)dd 2 ) (24) 
3W 

where cj) = r 2 /(Al 2 ) and j(4>) = y/cf> — r/(2l), with I an arbitrary length scale. Note that in the following <? M „ and 
associated quantities refer to a 2-dimensional metric. The action ([I]) becomes 

where 2G = G/l 2 , V{6) = 1/(2^), and 



2G J V R + "P J + 7 (25) 



1 = ~ = (26) 



Employing an ADM decomposition of the 2-dimensional metric in ( 24 ) 

ds 2 = e 2p (-a 2 dt 2 + (dx + Ndtf) (27) 
yields after some algebraic manipulation the Hamiltonian 

H = Jdx[^ (-M> + Anr 2 p M ) + (n- A (28) 

where M. is the Hamiltonian version of the Misncr-Sharpe mass function 

^ - e- 2 ? ((GII p) 2 - {4 >'f) + j -f (29) 
and we define the mass energy density pu via 

, 16' ( n 2 47^^ 2 e 2 ' , , A IGIV , , N 



The canonical momenta are given by 



n = ~(p'N + N'-p) (31) 



oG 

n p = (32) 
oG 



and we have defined 



47rr 2 

= C x ty-Ntf) (33) 



t = p'ii p -n' p + 0'n + n^' (34) 



It is important to note that the lapse, a, and shift, N, in the above parametrization are also the lapse and shift 
functions of the physical 4- metric as defined in ( 24 ) . 

Choosing the r as the spatial coordinate implies both the consistency condition = and 16' = j(6). Using this 
we obtain the shift in terms of the lapse: 
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Setting the diffeomorphism constraints strongly to zero, we eliminate 11^, yielding 



H = dx 



2p 



{-M 1 + Anr 2 p M ) + / dx(a^M)' 



J 



(36) 



for the partially reduced Hamiltonian upon including an appropriate boundary term. This partial coordinate choice 
allows many different fully gauge fixed theories, including Schwarzschild coordinates, Painleve-Gullstrand (P-G) co- 
ordinates and maximal slicings, among others. 

Specializing to P-G coordinates entails fixing g rr = 1, or e 2p = j ((/))■ The consistency condition p = fixes the 
lapse via 



(T 

a 



IGU^tp' 
2lGMj((f>) 



(37) 



determining p and a. Making use of the above conditions with eq (29) yields 

2GM _ (GTI P ) 2 



(38) 



with lip determined by the remaining Hamiltonian constraint. 

Finally, the field equations for the matter field independent variables ip and from the Hamiltonian ( 36 ) yields 
after some manipulation 



ip = a 



n, 



4Trr 2 £x V 



2GIM ,, 



[2gim 

\1W 



n, 



The quantities Ai and a as determined by the Hamiltonian constraint ( 38 ) 

- M' + 4:nr 2 p M - 



and the consistency condition ( 37 1 , giving 



with pm given by 



Airr p M = I 



a 
a 



n 2 



IGU^ 1 
2lGMj(<f>) 



(39) 
(40) 

(41) 
(42) 



l47rr 2 £y 



4?rr 2 £(X) 



'2GMI 

~1W 



(43) 



Recalling that 2G — G/l 2 , we note that the dynamical equations for the K-essence field and its conjugate momentum 
may be expressed in terms of four dimensional quantities only. 



As an example, for the choice ( 19 1 we have 



C 



1 



A" 



C 



1 



XX 



1 



2X 



-3/2 



In this case we can solve explicitly for X in terms of the phase space variables 



(44) 
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in P-G coordinates. 

For the purposes of our numerical simulations, we set the arbitrary length scale I to I = 1. Choosing the initial 
configuration for ip to be that given in equation (13 1, we first construct a fixed, spherically symmetric spatial grid 
in which to evolve the system, by specifying the number of spatial steps (we use 600), the minimum space between 
steps, and the maximum space between steps. Near the origin we require a finer grid. We therefore initialize the 
spatial grid with the minimum space between steps for a specified number of first steps nearest to the origin (we use 
100 steps). Once we are 100 steps from the origin, the spatial grid is gradually expanded so that the space between 
the grid points is at a maximum once we reach grid point 600. We select a minimum spacing of 1.0 x 10~ 3 and a 
maximum spacing of 1.0 x 10 -2 . After specifiying the initial configuration of the spatial grid, we then initialize the 
configuration of ip as specified above and we initialize the conjugate momentum Tl^p to be identically 0. 

After these initial conditions have been specified, we enter a loop, which proceeds at each timestep until a prede- 
termined and specfied maximal time is reached after which we anticipate no interesting evolution of the system. Note 
that while the spatial grid is constant at each timestep and only specified once, the timestep is not constant. It is 
determined during each execution of the loop. 

Inside the loop, we first calculate the necessary derivatives of tp and 11^, using finite difference methods. We then 
calculate C and its derivatives, as well as X and its first dervative in r. Next, we determine A4 and a on the spatial grid 
at the current timestep using a fourth order Runge-Kutta method. We then employ the condition that determines the 
location of the luminal apparent horizon and the sonic apparent horizon, as well as checking for Cauchy breakdown. 

At this point in the loop, data is written to file on selected executions of the loop. We then calculate the next 
timestep, according to the speed of a local ingoing geodesic. Finally, we evolve tp and 11^ in time using another fourth 
order Runge-Kutta method, and increase our total time by the calculated timestep. The loop then restarts. 

The loop continues to run until the total time has reached some predetermined and declared value as mentioned 
above. Alternatively, the program may cease to run due to a check performed on the values of C and its derivatives 
each time they are computed. If these values become divergent or imaginary, the program halts. This corresponds in 
the case of C\ to Cauchy breakdown, and in the case of £2 to the divergence of L^x ■ The program may also halt for 
£3 due to a check on the convergence of the Newton-Raphson solver. 
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